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Numerical test of Polyakov loop models in high temperature SU(2) 

Roberto Fiore a , Pietro Giudice a and Alessandro Papa a 
a Dipartimento di Fisica, Universita della Calabria 

& Istituto Nazionale di Fisica Nucleare, Gruppo Collegato di Cosenza, Italy 

We study the compatibility of effective mean-field models of the Polyakov loop for the deconfined phase of 
SU(N) pure gauge theories with lattice data obtained for the case of SU(2), in the temperature range T c ~ 4.8T C . 
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1. INTRODUCTION 

It has been suggested in several papers (see 
Ref. P) that the deconfined phase of SU(N) pure 
gauge theories could be described by an effective 
mean-field theory of the Polyakov loop, possess- 
ing global Z(N) invariance. Through this effective 
theory, a relation can be established between the 
pressure of the gluon gas and the Polyakov loop. 
If the phase transition is second order as for SU(2) 
or "weakly" first order as for SU(3), this effective 
theory can be written near the transition in terms 
of the first few powers of the Polyakov loop and 
of its complex conjugate. In this case, the rela- 
tion between pressure and Polyakov loop becomes 
very simple and its compatibility with lattice data 
can be easily tested. 

In this study, we have considered the case of 
SU(2) pure gauge theory on a 16 3 x 4 lattice with 
the standard Wilson action, in the temperature 
range T c ~ 4.80 T c . Although lattice effects are 
large for the Wilson action with N T = 4 sites in 
the time direction, the shape of the behavior of 
pressure and Polyakov loop with the temperature 
should not be different from the cases of larger 
values of N T , as seen in SU(3) 

2. LATTICE DETERMINATIONS 

The pressure of the gluon gas is given by 



P_ 



J_ 



/3 



= Nt 



ft 



d(3'[(S ) - (S T )} , (1) 



where Sq (St) is the action density at zero (non- 
zero) temperature, (3q is an arbitrarily chosen 
value, small enough that the integrand function 



at this point has become zero. Monte Carlo 
simulations were performed on 16 lattices for 
zero-temperature (typical statistics 30K), and on 
16 3 x 4 lattices for non-zero temperature (typical 
statistics 80K). Numerical results for N^[(Sq) — 
(St)] were interpolated by cubic splines before 
the numerical integration which led to the pres- 
sure (Fig. nj. As an estimate of the uncertainty 
for the pressure, we calculated also the integral 
by interpolating the data for N^[(So) — (St)} 
with the broken line connecting the la upper 
(lower) bound of each determination. The cor- 
respondence between f3 and the temperature has 
been established using the interpolating ansatz of 
Ref. .'! , which makes use of the known 0] critical 
couplings on lattices with A T =4, 5, 6, 8, 16. 

We considered both the charge- 1 and charge- 
2 Polyakov loops, given respectively by l\ = 
i(TrL(f)) and l 2 = ±(TrL(x) 2 > - [±(TrL(x))] 2 , 

with L(x) — rin T =i ^i(x, 714). We observe that l 2 
is Z(2)-invariant and is connected to the Polyakov 
loop in the adjoint color representation by Z a( jj — 
1 + 4Z2/3. In Fig. we show the behavior of if, 
If and I2 with (3. We observe that if goes lin- 
ear in the region 2.30 < f3 < 2.37 (corresponding to 
T c < T < 1.27 T c ) and in the region 2.45 < (3 < 2.80 
(corresponding to 1.65 T c < T < 4.80 T c ). More- 
over, l 2 goes to —3/4 in the confined phase, thus 
implying Z a( jj — > in that phase (for details on the 
behavior of Z a( }j across the transition, see Ref. [5])- 

3. POLYAKOV LOOP 
PURE GAUGE SU(2) 



MODELS IN 



Mean-field theory, dimensional analysis, Z(2) 
symmetry, reality of l\ in SU(2), power expan- 
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Figure 1. The three solid curves represent p/T 4 
and its uncertainty; the vertical lines represent 
the critical couplings on lattices with N T =A, 5, 6, 
8, 16 g|. 

sion in l\ imply the following simple form for the 
effective free energy: 



V 



^2 , 2 , °4 4 



if T 4 , b 2 > , 6 4 > . (2) 



The applicability domain of this model (called 
model A in the following) should be a region 
above T c , but not so close to T c that mean-field 
is spoiled, in which l\ is small enough to make 
if, If, ... terms negligible. The minimum of V is 
obtained for l\ = b 2 /b^ and leads to 

P_ _ Vmin _ h ^ 

According to this model, for constant 64, p/T 
should go linear with if. We find that the 
function (64/4) if fits the lattice data for the 
pressure in the region 2.33 < (3 < 2.37, i.e. 
1.11 T c <T< 1.27 T c , with 64 = 261.1(6.7) and 
x7(d.o.f.)=0.79 (see Fig. EJ. 

For high temperatures, one could expand the 
effective free energy in powers of (1 — If), thus 
getting 



VnT=[c-^q + ^it) t 



34 ,4 



which leads to 

P VjHT,min 



= c + b -±if 



(4) 



(5) 



Figure 2. if, If and l 2 vs (3 on a 16 3 x 4 lattice. 

There is compatibility of this functional form with 
lattice data for constant values of C and 64 in the 
region 2.60 < (3 < 2.80, i.e. 2.63 T C <T< 4.80 T c , 
with 64 = 19.7(4.8), C = 0.547(31) and 
X 2 /(d.o.f.)=0.18 (see Fig.|3J. 
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Figure 3. Comparison of the model A, for both 
low and high temperature regimes, with lattice 
data for the pressure. 

As a first variant of the model A, we consider 
the inclusion of the if term in the effective free 
energy (model B): 
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which leads to 
P_ _ Vmin = 



3 1 4 1 



(7) 



We find compatibility with the lattice data for 
the pressure over a wider region than in the case 
of model A, more precisely in the range 2.32 < 
(3 < 2.70, i.e. 1.07 T c < T < 3.56 T c , with 6 4 = 
350.3(6.1), b 6 = -1158(33) and x 2 /(d.o.f.)=0.73. 
A negative value for be would be problematic if 
the absolute value of l\ would be allowed to be- 
come arbitrarily large, which is not the case here. 
For high temperatures, using (1 — If) as expansion 
parameter we get 



P_ _ 

which agrees 
sure in the 
1.93T < T<4.80T, 



= C 



(8) 



with lattice data for the pres- 
region 2.50 < (3 < 2.80, i.e. 

ith 6 4 = 162(39), b 6 = 



-447(134), C = 0.258(65) and x7(d.o.f.)=0.15. 

Finally, we consider the model C obtained by 
model A with the inclusion of terms with the 
charge-2 Polyakov loop l 2 : 

leading to 



it+iilh) , (9) 
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For the high temperature version of this model, 
the only difference is an additive constant in the 
r.h.s. of the expression for p/T 4 . The compar- 
ison with lattice data shows that the inclusion 
of the terms with l 2 does not improve drasti- 
cally the quality of the fit in comparison with the 
model A. On the other side, the linear dependence 
of l 2 with If is roughly satisfied (x 2 /(d.o.f.);$2 
(see Fig. 0J) in both regions where the model 
A works, i.e. for 1.11 T c < T < 1.27 T c and for 
2.63 T C <T< 4.80 T c . This indicates that the be- 
havior of l 2 in T above T c is driven by the 
Polyakov loop l-y. The relatively large % 2 can be 
explained by the very small error bars both in l\ 
and in l 2 which make non-negligible higher pow- 
ers of l\ and I2 in the effective model. 
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Figure 4. l 2 vs l\ on a 16 3 x 4 lattice. 
There is a roughly linear dependence in two 
regimes: for 1.11 T c <T < 1.27 T c (corresponding 
to 0.038 < l\ < 0.061) and for 2.63 T C <T< 4.80 T c 
(corresponding to 0.135 < l\ < 0.182). 

4. CONCLUSIONS 

Lattice data show that p/T 4 has a roughly lin- 
ear behavior in a region centered around 1.2 T c 
and in a region centered around 3.5 T c ; in these 
regions also if exhibits a linear behavior, while 
li behaves linearly with l\. We have shown that 
both these evidences are in accord with simple 
mean-field effective models of the Polyakov loop. 
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